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Abstract 

The electromagnetic form factors of octet baryons are estimated within light cone 
QCD sum rules method, using the most general form of the interpolating current for 
baryons. A comparison of our predictions on the magnetic dipole and electric form 
factors with the results of other approaches is performed. 
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1 Introduction 



Electromagnetic form factors of nucleons provide information on their internal structures, 
i.e., about the spatial distribution of charge and magnetization of the nucleon. Nucleon 
electromagnetic form factors that are the functions of only four-momentum transfer squared 
Q 2 are described by Dirac F\(Q 2 ) and Pauli i^Q 2 ) form factors which are related to the 
electric and magnetic dipole form factors Ge{Q 2 ) and Gm(Q 2 ) as, 



Obviously, in the limit Q 2 — > the form factors Ge and Gm correspond to the charge 
and magnetic moments of the nucleon, while F\ and F 2 describe the charge and anomalous 
magnetic moments of the nucleon. 

The study of electromagnetic form factors of hadrons receives a lot of attention during 
the past decade. Recent experiments on the nucleon form factors using the polarized elec- 
tron beam and polarized protons, which are presented in detail in [1], allow more accurate 
measurements of the nucleon form factors at higher values of the momentum transfer. In 
the polarization measurements it is observed that the ratio G^(Q 2 )/G^(Q 2 ) can not be 
determined by the simple dipole form G D (Q 2 ) = (1 + Q 2 /Ql)~ 2 with Ql = (0.71 GeV) 2 
[2-4]. The neutron form factors that are measured up to Q 2 = 3.4 GeV 2 recently can 
provide detailed comparison of the proton and neutron form factors [5,6]. 

Considerable progress has also been achieved on the electromagnetic excitation of nu- 
cleon resonances during last years. The cross sections and photon asymmetries for the 
photo production of the pion and rj mesons are measured at MAMI at Mainz, ELSA at 
Bonn, LEGS at Brookhaven, and GRAAL at Grenoble. Moreover, a large amount of data 
has already been collected for the A(1232) excitation and single Q 2 data points are obtained 
for the longitudinal and transversal form factors of the p — > A(1232), Pn(1440), S f n(1535), 
_Di3(1520), etc., whose results are all given in [7]. These progresses in experiments open a 
way to real possibility of measuring the electromagnetic form factors of the octet baryons 
in near future. 

In the present work we calculate the electromagnetic form factors of the octet baryons 
within the light cone QCD sum rules (LCSR) method by employing the general form of the 
interpolating currents. It should be noted that this problem is studied in the same method 
for the Ioffe current alone in [8-10]. It should also be reminded to the interested reader that 
the nucleon electromagnetic form factors are calculated for the Ioffe and general currents 
in [11] and [12], respectively. 

The plan of this work is as follows. In section 2 we introduce the correlation function 
which we shall use in calculating the electromagnetic form factors of the octet baryons, and 
discuss how the interpolating currents of the octet baryons are related to each other. In 
section 3, the light cone QCD sum rules for the electromagnetic form factors are obtained in 
the case when the correlation functions are calculated in terms of the main nonperturbative 
input parameters, namely in terms of distribution amplitudes (DAs) of the octet baryons. 
The last section contains the details of the numerical calculations of the electromagnetic 
form factors of the octet baryons. 




F 2 (Q 2 ) , 



G M = F l (Q 2 ) + F 2 (Q 2 ) . 



(1) 
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2 LCSR for the electromagnetic form factors of octet 
baryons 

In order to obtain the LCSR for the electromagnetic form factors of the octet baryons we 
start by considering the following vacuum-to-one-octet baryon correlation function, 

U,(p,q) = i J d'xe^ (0 \T{ V (0)fJ(x)}\B(p)) , (2) 

where r] is the interpolating current of the octet baryon, jf{x) is the electromagnetic current, 
ji is the vector Lorentz index, T is the time ordering operation, and B{p) is the one particle 
baryon state with momentum p. 

The most general forms of the interpolating currents for the octet baryons can be written 

as, 

= V2e abc J2 [(u aT CA[s b y e 2 d c + (d aT C A[s b ) A l 2 u c 
i=i 

rf + =2e abc (u aT CA[s b ^Aiu c , 

rf~ = 2e abc [d aT CA[s b ^A e 2 d c , 
-° s+ / \ 

T] =7] [U -H- S) , 

rf~ =rf-(d<r>s) , (3) 

where A\ = I, A\ = A\ = 75, A\ = (3. 

The interpolating current of the A baryon can also be obtained from that of S° baryon 
in the following way [13]: 

2r]^\d <H> s) + rf Q = -VSr) A , or, 

2rf\u^ s)-rf° = -V3 VA . (4) 

Our primary aim is the calculation of the phenomenological part of the correlation 
function (2). According to the standard procedure, in order to obtain the physical part of 
the correlation function of the octet baryons we insert a full set of baryons into Eq. (2). 
Separating the contribution of the ground state baryon we get, 

(0 \rj\ B(p - q ))(B(p-q) B{p)) 
m 2 B -(p- qf 

where • • • represents the contributions of the higher states and continuum. 

The matrix element (0 \q\ B{p — q)) appearing in Eq. (5) are determined as, 

(0|i/|B(p-?))=Aflii(p-g) , (6) 

where As is the residue of the members of the octet baryons. The hadronic matrix element 
with the electromagnetic current is determined in terms of three independent form factors 
Fi, F 2 and F 3 in the following way, 

(B(p-q)\jf!\B(p)) =u(p-q)[F 1 (q 2 ) 1 ,-t^F 2 (q 2 ) + ^F 3 (q 2 )]u(p) . (7) 



n„(p,g)= v S 1* +••• > ( 5 ) 
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From conservation of the electromagnetic current we get F 3 (q 2 ) = 0. Taking Dirac equation 
into account, one can show that the general decomposition of the correlation function (2) 
contains six independent amplitudes in the presence of the electromagnetic current, 



n M (p, q) = + u 2 pj + n 3 7/, + n 4 7^ + n 5 g M + n 6 g^ u(p) . 



(8) 



Using the definitions given by Eqs. (6) and (7), we get the following expression for the 
hadronic part, 



n M (p,g) 



A 



B 



+ 



m%-{p- q) 
- 2F 1 (q 2 ) - F 2 (q 2 ) 



2F 1 (q 2 )p, + ^-pj + fi^g 2 ) + F 2 (q 2 )}^ 



F 2 (q 



m B 

2\ 



2m i 



(9) 



Equating the coefficients of each Lorentz structure in Eqs. (8) and (9) we get the sum rules 
for the form factors. In order to perform numerical analysis we need expressions of the 
invariant functions Ilj, (i — 1, • • • , 6) from QCD side. 

The calculation of the invariant functions Ilj from QCD side is carried out when the 
external momenta p — q and q are taken in deep Eucledian space, i.e., (p - g) 2 < and 
q 2 <C 0, which is necessary to perform operator product expansion (OPE) near the light 
cone x 2 ~ 0. The OPE result can be obtained as the sum over octet baryon distribution 
amplitudes (DAs) of growing twist, which are the main non-perturbative inputs of the 
LCSR. 

As has already been noted, the DAs of S, S and A are investigated in [8-10]. The DAs of 
the octet baryons with J p = l + /2 are defined from the matrix element of the three-quark 
operator between the vacuum and the baryon state \B(p)), whose form is given as, 



(0 \qi a (aix)q b 2fS (a 2 x)q c 3l (a 3 x) \ B(p)} , 



(10) 



where a, (3, 7 are the Dirac indices, a, b, c are the color indices, and are positive numbers 
satisfying a± + a 2 + 03 = 1. Using the Lorentz covariance, as well as spin and parity of the 
baryons under consideration, the matrix element (10) can be decomposed as, 

Ae abc (0 \ q a la ( ai x)q b 2fi (a 2 x) q c 37 (a 3 x)\ B{p)) = £ T^t (^(p)) , 



(11) 

where are certain Dirac matrices, Ti = Si, Vi, Ai, Vi and % are the DAs which do not 
have definite twists. The DAs with definite twists are determined from, 

4e afec (0 ^U^qwi^qt^a^ B{p)) = £ Ftff (^(p)) , (12) 

i 

where Fj = Si,Pi,Ai,Vi, and T { . The Relations among these two sets of DAs are given 
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as, 



5i = 


= 5! , 




(2P 


x) 


s 2 = 


S\ — s 2 , 


Pi = 


-Pi , 




(2P 


x) 


v 2 = 


p 2 - Pi , 


Vi = 


-■v 1 , 




(2P 


x) 


v 2 = 


v 1 -v 2 -v 3 , 


2V 3 


= ^3, 




(4P 


x) 


V 4 = 


-2V 1 + V 3 + V 4 + 2V 5 , 


(4P 


x)V 5 = 


v 4 -v 3 , 


(2P 


x) 


2 v 6 = 


= -Vx + V 2 + V 3 + V 4 + V 5 - V 6 , 


A = 


-A 1 , 




(2P 


x) 


A 2 = 


-A l + A 2 - A 3 , 








(4P 


x) 


A 4 = 


-2A 1 -A 3 -A 4 + 2A 5 , 


(4P 


x) A 5 = 


A 3 - A 4 , 


(2P 


x) 




= A 1 - A 2 + A 3 + A 4 - A 5 + A 6 , 


Ti = 


Ti , 




(2P 


x) 


r 2 = 


Ti + T 2 - 2T 3 , 


2T 3 


= T 7 , 




(2P 


x) 


% = 


Ti - T 2 - 2T 7 , 


(2P 


x)% = 


_ Tl + T 5 + 2T 8 , 


(2P 


x) 


2 % = 


= 2T 2 - 2T 3 - 2T 4 + 2T 5 + 2T 7 + 2T f 


(4P 


x)T 7 = 


T 7 -T 8 , 


(2P 


x) 


2 T 8 = 


= -7\ + T 2 + T 5 - T 6 + 2T 7 + 2T 8 . 



The complete decomposition of the DAs in Eq. (11) in terms of Si, Vi, Ai, Vi and % func- 
tions, as well as the explicit expressions of DAs, can all be found in [8-1 1] . 

Omitting the details of calculations of the theoretical part and equating the coefficients 
of the Lorentz structures p^, p^ from hadronic and theoretical parts, and performing the 
Borel transformation and continuum subtraction over the variable (p — q) 2 , we get the 
following sum rules for the form factors, 



Fi{Q 2 ) 



2\ 



p 6 (x) 



X 



X 



2 M 2 2x 3 M* J 



r e x 



2 xm\. 

tH V- 

M 2 M 1 



+ 



1 

+ 7^" 



Pe(^o) 



B y. Jx 

Q 2 + x\m\ 

2 

2 (Q 2 + x 2 m 2 B ) \dx x (Q 2 + x 2 m 2 B )M 2 



2x (Q 2 
d 



x 2 m 2 B )M 2 
Pe(^o) 



-( S0 -m|)/M 2 



(13) 



F 2 (Q 2 ) = 2m B F 1 (Q 2 ) {p 2 {x) ->• p' 2 (x), p 4 (x) -)• p' 4 (x), p 6 (x) ->• , (14) 

x 

where M is the Borel parameter, s = — Q 2 + (1 — x)m 2 B , x is the solution of the equality 

x 

s = So, rriB is the mass of the members of the octet baryons and x — 1 — x. The factor L 
in Eq. (13) is the normalization factor whose value for the members of the octet baryons 
is determined as, 



f 1 
2 



V2 



for E+E", 

for S°, 
for A . 



i — i , i — i , 



(15) 



The explicit expressions of pi and p\ for S + , S° and A baryons are presented in the Ap- 
pendix. 
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3 Numerical analysis 



As has already been mentioned, the main nonperturbative inputs of LCSR are the baryon 
DAs. Here we would like to make the following remark about the expressions of the DAs for 
the A, E and S baryons. In [14], the DAs for nucleons were extended up to next-to leading 
order in conformal spin and the expressions of the nucleon DAs of twist-3 up to next-to-next 
to leading conformal spin is found in [15]. As a result of these two works it is obtained that 
the nucleon form factors are sensitive to the higher conformal spin contributions. For other 
members of the octet baryons similar calculations are not yet done and deserves a detailed 
study. In present work, we consider the DAs for the A, E and S baryons without these 
contributions, whose expressions can be found in [8-11]. The parameters appearing in the 
expressions of the DAs are estimated from the analysis of the sum rules given in [10-12]: 



fs 


= (9.9 ±0.4) x 10~ 3 


GeV 2 , 


Ai 


= -(2.1 ±0.1) x 10" 


2 GeV 2 


A 2 


= (5.2 ±0.2) x 10~ 2 


GeV 2 , 


A 3 


= (1.7 ±0.1) x 10~ 2 


GeV 2 , 


h 


= (9.4 ±0.4) x 10~ 3 


GeV 2 , 


Ai 


= -(2.5 ±0.1) x 10" 


2 GeV 2 


A 2 


= (4.4 ±0.1) x 10~ 2 


GeV 2 , 


A 3 


= (2.0 ±0.1) x 10~ 2 


GeV 2 , 


h 


= (6.0 ±0.3) x 10~ 3 


GeV 2 , 


Ai 


= (1.0 ±0.3) x 10~ 2 


GeV 2 , 


A 2 | 


= (0.83 ±0.05) x 10 


2 GeV 2 


A 3 | 


= (0.83 ±0.05) x 10 


2 GeV 2 



The remaining input parameters of the LCSR are the continuum threshold s , the 
Borel parameter M 2 and the auxiliary parameter j3 that appears in the expressions of the 
interpolating currents of the octet baryons. 

In our numerical calculations we use the values sq = 2.5 GeV 2 , sq = 3.0 GeV 2 and 
So = 3.2 GeV 2 for the continuum threshold, obtained from mass sum rules analysis in [16], 
for the A, E and S baryons, respectively. 

The Borel mass parameter M 2 is another auxiliary parameter of the sum rules. There- 
fore the "working region" of M 2 should be found, where the form factors are practically 
independent of it. The lower limit of M 2 can be obtained by requiring that the higher 
states and continuum contributions to the sum rules constitute, maximally, about 40% of 
the total result. The upper limit of M 2 can be determined by demanding that the operator 
product expansion should be convergent. Our calculations show that the region in which 
the aforementioned two conditions are properly satisfied, are: 1.2 GeV 2 < M 2 < 3.0 GeV 2 
for E and A baryons; and 1.4 GeV 2 < M 2 < 3.5 GeV 2 for S baryons. In further numerical 
analysis, we use M 2 = 2.0 GeV 2 for E, S and A baryons. 
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The residues of the octet baryons are calculated in [16] and we shall use these results in 
our numerical analysis. Furthermore, it should be noted that, from experimental point of 
view, it is more convenient to study the Sachs form factors Ge and Gm that are given in 
Eq. (1). 

The Q 2 dependence of the magnetic and electric form factors for S + , S and A, baryons 
are shown in Figs. (l)-(6). In order to get "good" convergence of the light cone expan- 
sion and reliable results from the LCSR, sufficiently large Q 2 is needed. In our numerical 
calculations we consider the lower limit of Q 2 as Q 2 = 1 GeV 2 , where above this point 
the higher twist contributions are suppressed. On the other side, the higher resonance and 
continuum contributions become small enough when Q 2 < 8 GeV 2 . For this reason, we 
perform numerical analysis in the region 1 GeV 2 < Q 2 < 8 GeV 2 . 

In odd-numbered Figs.: (1), (3), (5) (even-numbered Figs.: (2), (4), (6)) we present 
the dependence of the magnetic dipole form factors Gm(Q 2 ) (electric form factor Ge{Q 2 )) 
on Q 2 , at fixed values of s and M 2 chosen from their working regions, and at several fixed 
values of the arbitrary parameter (3, for S + , H~ and A baryons, respectively. 

• In the case of S + , these form factors get positive (negative) values for negative (pos- 
itive) values of the parameter (3. 

• The situation for S~ is contrary to the behaviors of the form factors of the S + , i.e., 
the values of Gm and Ge are positive (negative) when the parameter j3 is positive 
(negative). 

• In the case of S° baryon, the form factors exhibit similar behaviors as the form factors 
of S + baryon. 

• It is observed that the form factor Ge of the S° baryon changes its sign practically 
at all considered values of (3. The zero values of Ge depend on the values of the 
arbitrary parameter /3. But the values of Ge are quite small, whose maximum value 
is about G E (Q 2 ) = 0.06. 

• It is interesting to observe that at j3 — 3 and (3 = 5, Gm for the S° changes its sign, 
while for the other values of (3 it is always negative. 

• For positive (negative) values of [3 the magnetic dipole form factor Gm{Q 2 ) for A 
baryon attains at positive (negative) values. 

• The situation for electric form factor Ge{Q 2 ), however, is slightly different. Namely, in 
the case of Ioffe current for which (3 = —1, Ge(Q 2 ) becomes negative whose magnitude 
is negligibly small. 

• When the form factors of the S~ baryon are considered we see that Gm changes its 
sign only at /3 — 5, while for all other values of (3 it gets at only negative values. On 
the other hand the form factor Ge is positive (negative) for all positive and negative 
values of (3. 

We now compare our results on Q 2 dependence of the magnetic and charge form factors 
with the ones existing in the literature. These form factors are discussed within the LCSR 



6 



method for the Ioffe current (/3 = — 1) in [9, 10], within the framework of the fully relativistic 
constituent quark model in [17], in the framework of the covariant spectator quark model 
[18] and lattice QCD [19]. When we compare our predictions on the form factors with the 
results of the above-mentioned works we obtain that, at j3 — — 1, our predictions on are 
very close to the results of [9, 10], and [17], except for the S° baryon. Our predictions for 
the magnetic form factor Gm agree within the errors with the existing results. The small 
differences among the predictions can be attributed to the different values of the input 
parameters used in the numerical analysis. 

As has already been noted, the interpolating currents of the octet baryons contain also 
the auxiliary arbitrary parameter j3. Obviously, the physically measurable quantities should 
be independent of this parameter. In order to find the working region of the parameter (3 
we demand that the form factors are practically independent of it. As an example, in Figs. 
(7) and (8) we present the dependence of Gm(Q 2 ) and Ge(Q 2 ) on cos 9 (j3 = tan 9) for S + 
baryon at fixed values of s = 3.0 GeV 2 and M 2 = 2.0 GeV 2 , for two fixed choices of Q 2 , 
namely, Q 2 = 2.0 GeV 2 and Q 2 = 3.0 GeV 2 . We see from these figures that, in the region 
—0.3 < cos# < 0.3, Ge and Gm show very weak dependence on (3. In other words, the 
working region of the parameter (3 for the E + baryon is —0.3 < cos# < 0.3. 

We perform similar analysis for all other members of the octet baryons and find out 
that the region —0.2 < cos6> < 0.2 is the common working region to them as well. It 
should be noted here that Gm(Q 2 ) for S° and S~ baryons exhibit stability in the range 
—0.2 < cos# < 0.2. Also note that /3 = — 1 point, which is the Ioffe current corresponding 
to cos# = —0.71, belongs to the region where the predictions for the form factors are not 
reliable. Choosing the values of M 2 and f3 from the relevant working regions, and from a 
comparison of our predictions on the form factors with the results of the above-mentioned 
works we see that, 

• Predictions of all works for Gm{Q 2 ) are very close to each other within the error 
limits. 

• Our predictions on Ge{Q 2 ) agree with the results of other approaches, except for the 
S° and A baryons. In these cases our results are very close to the predictions of the 
lattice QCD, while considerable disagreements are observed with those obtained in 
[17] and [18]. 

The results obtained in this work can be improved by taking into account the 0(a s ) 
corrections to the distribution amplitudes, and more accurate values of the input parameters 
entering the sum rules. 

In conclusion, in the present work we have studied the charge Ge{Q 2 ) and magnetic 
dipole Gm{Q 2 ) form factors within the LCSR method by using the most general for of the 
interpolating currents for the octet baryons. We have compared our predictions on these 
form factors with the results existing in literature that were obtained in framework of the 
relativistic quark model, lattice QCD and LCSR for the Ioffe current. 
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Appendix 

In thhis Appendix we present the explicit expressions of the functions p 2 {p 2 ), Pa{p'a) an d 
Pg(p'o) entering to the sum rules for the form factors Fi(Q 2 ) and F 2 (Q 2 ), 



ptW 



+ 



+ 
+ 
+ 



4e u m|+(l + f3)x(m\ + x 2 + Q 2 )B 6 (x) 

4e s m|+|m|+m s (l - /3)x 2 C 6 + (1 + (3) m s +(m|+x 2 + Q 2 )xB 6 

- m s (Q 2 5 6 + 2m| + x 2 5 8 )]}(x) , 

pf{x) = e u m s+ { - 2m|+a;[(l - f3)(d 6 +D 6 ) - (1 + p)(2B 6 - 3£ 8 )] (x) 
+ (1 - p) [m|+x 2 (£>4 - 3D 5 - C 4 + 3C 5 ) + 2Q 2 (D 2 + C 2 )l (x) 
+ (1 + p) \Q\B 2 + 55 4 ) - ml + x 2 (2H l - 2E X + B 2 -B 4 + QB 5 + 12£ 7 )] (x) 

- 2m| + x ^ dx 3 [(1 - /?)« - V X M ) + 3(1 + /3)7f] (x, 1 - x - x 3 , x 3 ) } 

e d m s+ 1 - 2m| + x [(1 - p)(C 6 -D 6 ) + 2(1 + p) S 8 ] (x) 
(1 - 0) [ - m|+x 2 (5 4 - 5 5 + C 4 - C 5 )] (x) 
(1 + p) [2Q 2 {B 2 + B 4 ) - 4m| + x 2 (5 5 + 2S 7 )] (x) 

2m| + x jf* tfax [(1 - /3)(Af + Vf) + 2(1 + /3)Tf] (x 4 ,x, 1 - x x - x)} 

e s m E +|2m s +(l + /3) m Ti +x(2B & -B 8 ) - m s B & (x) 
(1 - /3) [2(m 2 +x 2 C 5 + Q 2 d 2 ) - m s+ m s x(2C 2 - C 4 - C 5 )l (a:) 
- (1 + (3) Q 2 (B 2 - 3B 4 ) + ml+x 2 (B 2 -B A + 2B 5 + 45 7 ) - 4m s+ m s x(5 4 - S 5 )] (x) 

- 2m|+(l + (3)x J dxi T^(xi, 1 — x\ — x, x)| , 

pf(x) = 2e„m s+ x / cb 3 [(1 - p)(A x + 2A 3 - K + 2V 3 ) 

- (1 + /3)(Pi + 5i + 3Ti - 6T 3 )] (x, 1 - x - x 3 , x 3 ) 

+ 2e d m s+ { [(1 - /3)(5 2 - C 2 ) + (1 + /3)(5 2 - £ 4 )] (x) 

- x J dx x [(1 - p){A x +A 3 + V 1 - V 3 ) + 2(1 + /3)(Tx - 2T 3 )] (x 1; x, 1 - x x - x)} 
+ 2e s {m s+ [(1 - /3)C 2 + (1 + /3)(B 2 - 5 4 )] (x) 



+ 
+ 
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+ f d Xl \(1 - f3)(mx+xV 3 + rrisVi) 
Jo 1 

- (l + /3)(m E+ x(Pi + S , i + T 1 -2T 3 ) + m s T 1 )](xi,l-Xi-x,x)} , 

P? + {x) = -4e tt m|+(l + /3)(m|+x 2 + Q 2 )5 6 (x) 

- 4e s m|+{m s+ m s (l -/3)x(7 6 + (1 + /3) (m|+x 2 + Q 2 )B & + m s +m s x(-B 6 — 2_B 8 )j j(x) , 

P? + (x) = e u m|+{ -3(l + /3)B 6 (x) 

+ x[(l - /3)(2 J D 2 - D 4 + 3,D 5 + 2(7 2 + C 4 - 3C 5 ) 
+ 2(1 + /?)(#! - + B 2 + 25 4 + 35 5 + 6S 7 )] (x)} 
+ e,m| + {2(l + /3)5 6 (x) 



+ x 



(1 - p)(D 4 -D 5 + C 4 - C 5 ) + 2(1 + (3){B 2 + B 4 + 2B 5 + AB 7 ) 

+ 2(1 - 0) jf « + V^)^, X, 1 - X! - x)} 

{ -5m s+ (l + /3)5 6 (x) 
- 2 (1 - /3) (m s+ xC 5 - (m s+ x + m s )C 2 ) - (1 + £) (m E+ x(5 4 + 5 5 + 25 7 ) - m s (B 2 + B 4 ) 

+ 2m E+ (l-£) jf rfxi(Af -\/ 1 M )(x 1 ,l-x 1 -x,x)} , 

' 2 s+ (x) = 2e d (l - /3) / {A 1 + ^(x^x, 1 - xx - x) 
Jo 

-2e 8 (l-/3) dx 1 Vi(x, 1 - x - x 3 ,x 3 ) , 
./o 

pf (x) = 4e u m| (l + /3)x(m| x 2 + Q 2 )B 6 (x) + 4e d m| (l + /3)x(m| x 2 + Q 2 )B 6 (x) 

(1 - /3)x 2 C 6 + (1 + (5) U s ox(m 2 x 2 + Q 2 ) B 6 



(x) 



Pi 



m s (Q 2 B 6 + 2m| x 2 5 8 )]}(x) , 

) = e u m S o{ -2m| x[2(l-/3)(7 6 -(l + /3)(25 6 -55 8 )](x) 
+ [2(1 - (m| x 2 (£> 5 - + 2C 5 ) - Q 2 (D 2 - C 2 )) 

+ (1 + /3) (Q 2 (35 2 + 75 4 ) + ml x 2 (2H 1 - 1E X -B 2 + B 4 - 10B 5 - 20£ 7 ))] (x) 
- 2m| x J dx 3 [2(1 - f3)V™ + 5(1 + /3)Tf] (x, 1 - x - x 3 , x 3 )} 
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+ e d m s o| - 2m| x[2(l - p) C 6 - (1 + /3)(2B 6 - 5B 8 )] (x) 
+ [(1 - P) ( - 2m| x 2 ( J D 5 + C 4 - 2C 5 ) + g 2 (5 2 + C 2 )) 

+ (1 + /3)(g 2 (3fi 2 + 75 4 ) - mloX 2 (2H l - 2E X + B 2 -B 4 + 10B 5 + 20£ 7 ))] (x) 
-2ml xJ^ dx x 2(1-/3)^ + 5(1 + /?)^ (x 1; x, 1 - x x - x)} 
+ 2e s m s o|2m s o(l + (3) m^ox(2B 6 -B 8 ) - m s B 6 (x) 



+ 



(1 - (3) (2(mloX 2 C 5 + Q 2 C 2 ) - m^m s x{2C 2 -C 4 - C 5 )) 



- (1 + 0)1 Q\B 2 - W 4 ) + mi z0 x\B 2 -B 4 + 2B 5 + 4E 7 ) - 4m E om s x(B 4 - 5 5 



pf\x 



2m 2 z0 {l + P)x J dxi T™(xi, 1 — x± — x, x)j , 

) = -2e u m s o{ [(1 - £)(£> 2 + C 2 ) - (1 + (3)(B 2 - B 4 )j (x) 
+ x dx 3 [(1 - P)(A 3 + 2V! - 3V 3 ) - (1 + P)(Pi + Si - 57\ + 10T 3 )] (x, 1 - x - x 3 , a; 3 )} 
+ 2e d m s o{ [(1 - /3)(5 2 - C 2 ) + (1 + /3)(B 2 - 5 4 )] (x) 

+ rfxi [(1 - /3)(A 3 - 2Vi + 3V 3 ) - (1 + /3)(P 4 + 5i + 5Tx - 10T 3 )] (xi,x, 1 - x 4 - x)} 

+ 4e s {m s o [(1 - /3)C 2 - (1 + /3)(B 2 - B 4 )j (x) 

+ y dxi |(1 - 0)(m^oxV 3 + m s Vi) + (1 + f3) 2m s oxT 3 - (m s ox + 2m s )TiJ |(xi, 1 - Xi - x, a 

P? ^) = -4e u m| (l + /3)(m| x 2 + Q 2 )£ 6 (x) - 4e d m| (l + /3)(m| x 2 + Q 2 )5 6 (x) 

- 8e s m| |m s om s (l - 0)xC & + (1 + /?) (m| x 2 + Q 2 )B e + m E om s i(B 6 - 25 8 ) |(x) , 

P 4 E °(z) = -e„m| {(l +^)i 6 (x) 

+ 2x [(1 - /3)(I> 2 + A, - C 2 - C 4 + 2C 5 ) + (1 + /3)(#i - E ± - 2B 2 - 3B 4 - 5B 5 - 10S 7 )] (x) 

+ 2(1 - 0) ^ rfx 3 « - V; M )(x, 1 - x - x 3 , x 3 ) } 

+ e(i m 2 { - (l + /3)5 6 (x) 

+ 2x [(1 - f3)(D 2 + D 5 + C 2 + C 4 - 2C 5 ) + (1 + /?)(#! - E x + 25 2 + 3£ 4 + 5B 5 + 105 7 )] (x) 
+ 2(1 -P) j d Xl {£f + y i M )(x 1 ,x,l-x 1 -x)} 

- 2e s m s o i 5m s o (1 + /3)B 6 (x) 
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(1 - 0) (mv>xC 5 - (m s ox + m s )C 2 ) -(1 + 0) (m^x(B 4 + B 5 + 2B 7 ) - m s (B 2 + £? 4 )) (x) 
+ 2m^a(l-0) J drn^^i,! - xi , 



p 2 E °(:r) = -2e u (l - / (fa 3 (A 1 - V 1 )(x, 1- x - x 3 , x 3 ) 

Jo 

rx 

+ 2e d (l-0) / dx 1 (A 1 + V 1 )(x 1 ,x,l-x 1 -x) 
Jo 

rx 

-Ae s (l-0) dxi V 4 (x, 1 - x - x 3 ,x 3 ) , 
Jo 

p$(x) = -12e u m A (l + /3)x(mfc 2 + Q 2 )£ 6 (:r) - 20e d m|(l + /3)x(m 2 v x 2 + Q 2 )5 6 (:r) 
+ 8e s m 2 A ^2m s Q 2 B 6 - m 2 A m s x 2 (C 6 - 2B 8 ) - m A (l + 0)x 2(m\x 2 + Q 2 ) B 6 

+ m A m s xD 6 ^ , 
fi(x) = 3e u m A {2m 2 A x [2(1-/3) d 6 - (1 + 0){2B 6 - bB^x) 



+ 
+ 



2(1 - 0) [m 2 A x 2 (D 5 -C 4 + 2C 5 ) - Q 2 (D 2 - C 2 
+ (1 + 0) (Q 2 (3B 2 + 7B 4 ) + m 2 A x 2 (2H 1 - 2E X - B 2 + B 4 - 10B 5 - 20S 7 ))] (x) 

+ 2m A x jf cfe 3 [2(1 - 0)V? + 5(1 + ^)T 1 M ] (x, 1 - x - x 3 , x 3 )} 

+ e,m A |2m 2 :r [2(1 -p)C 6 -(l + 0)(1OB 6 - 9B S ) 
2(1 - 0) (m 2 A x 2 (D 5 + C 4 - QC 5 ) - Q 2 (D 2 + 5C 2 ) 

(1 + /3) (Q 2 (B 2 - 19B 4 ) + m\x 2 (2H l - 2E X + 55 2 - 5B 4 + 185 5 + 365 7 ) 

+ 2m l x J q dx i 2 (! - /W M + 9(! + £) T i M (xi,x, 1 - x i ~ x )} 
+ 2e s m A \Am A m A (l - 0)xC 6 - m A (l + 0)x(2B Q - 3B$) + m s B^ (x) 
+ [2(1 - 0) (m 2 A x 2 (D 5 + C 4 - 3C 5 ) - Q 2 (D 2 + 2C 2 )) 
+ 2m A m s :r(2C 2 -C A -C 5 - AB 4 + 4B 5 ) 

- (1 + 0) (2Q 2 (B 2 + 5B 4 ) - mlx 2 ^ -E ± + B 2 -B 4 + 6B 5 + 12B 7 ) 

- m A m s x(2D 2 + D 4 + l) 5 )) (x) 

+ 4m 2 A xJ^ dx! (1- 0)Vf* + 3(l + 0)Tf* (x u 1 - Xi - x, x)} 
p£(x) = 6e u m A { [(1 - /3)(L> 2 + <7 2 ) - (1 + /3)(5 2 - S 4 )] (x) 
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+ x j dx 3 [(1 - P)(A 3 + 2V l - 3V 3 ) - (1 + /3)(Pi + S-x - 57\ + 10T 3 )] (x, 1 - x - x 3 ,x 3 )} 

- 2e d m A { [(1 - 0)(D 2 + 3C 2 ) - 3(1 + (3)(B 2 - S 4 )] (x) 

+ x^ rfx! [(1 - /3)(A 3 - 2Vi + 7F 3 ) - (1 + P)(Pi + Si + m - 18T 3 )] (x 1; x, 1 - Xl - x)} 

- 4e s {m A (l-/3)5 2 (x) 

+ J dx x [m A (l - /3)x(A 3 - 2Vx + 4\/ 3 ) - (1 + 0) {m K x{P 1 + S 1 + 67\ - 12T 3 ) + m s A^j 
+ 2m s (V 1 -2T 1 ) (xi, 1 - xi - x,x)| , 



p 6 A (x) = 12e u m A (l + f3)(m 2 A x 2 + Q 2 )B 6 (x) + 20e d m A (l + f3)(m 2 A x 2 + Q 2 )B 6 (x) 

+ 8e s m A {m A m s x(C 6 +5 6 - 2§ 8 ) + (1 + 0) 2(m A x 2 + Q 2 ) B 6 + m A m s xD 6 |(x) , 

pl A (x) = 3e u m A {(l + /3)il 6 (x) 

+ 2x [(1 - /3)(L> 2 + A, - C 2 - C 4 + 2C 5 ) 

+ (1 + £)(#i -E x - 2B 2 - 3B A - 5B 5 - IOS7)] (x) 

+ 2(1 - 0) jf ^ 3 (Af - V; M )(x, 1 - x - x 3 , x 3 )} 

+ e,m A {21(l + /3)5 6 (x) 

- 2x[(l - p)(D 2 + D 5 + 5C 2 + C 4 - 6C 5 ) 

+ (1 + P)(H l -E ± + 2B 2 + 7B 4 + 9B 5 + I8S7)] (x) 

- 2(1 -P)J dx l { A l ~ 3V™){x 1 ,X, 1 - X! - x)} 

+ 4e s m A |3m A (l + /3)5 6 (x) 

- (l-/3)(m A x(D 2 + D 5 + 2C 2 + C 4 -3C 5 )^ 

+ (1 + 0) [m A x(H l -E 1 + 2B 2 + 4B 4 + QB 5 + 12B 7 ) + m s 5 2 ) + 2m s (C 2 - B 2 - B 4 ) 

- m A (l — j3) J dxi A^(xi, 1 — xi — x, x)| , 

p 2 A (x) = 6e u (l - /3) / rfx 3 (Ai - V 1 )(x, 1 - x - x 3 ,x 3 ) 



-2e d (l-/3) f dx 1 (A 1 -3V 1 )(x 1 ,x,l-x 1 -x) 
Jo 
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4e s (l-f3) dx 1 A 1 (x,l-x-x 3 ,x 3 ), 
Jo 



In the above expressions for pi and p[, the functions J-(xi) are defined in the following 
way: 

dx[ J dx 3 F(x' 1: 1 - x\ - x 3 , x 3 ) , 

JF{xi) = / da^ / dx 1 j dx 3 J r (x 1 , 1 — x 1 — x 3 ,x 3 ) , 
ii </i Jo 



/x 2 ^ r 1 -~ x 2 
dx 2 / dxi J-(xi, x 2 , 1 — x\ — x 2 ) , 

f* X 2 ^ 2 /* 1 "^2 



^{x 2 ) = / / ^2 / .F(Xi, X 2 , 1 — Xi — X 2 ) ' 

il ii Jo 
^ /" X3 / /" 1_:E » 

•T""^) = / dx 3 / dxi J-(x\, 1 — xi — x 3 , x 3 ) , 

•F(x 3 ) = / dx' 3 / cfo 3 ' / dx 1 F(x u l - x x - x" 3 ,x" 3 ) , 

Jl Jl Jo 

and the definitions of the functions B iy Cj, A, £1 and /Ti are given as 

5 2 = Ti + T 2 - 2T 3 , 

5 4 = Ti - T 2 - 2T 7 , 

£? 5 = -7\+T 5 + 2T 8 , 

5 6 = 27\ - 2T 3 - 2T 4 + 2T 5 + 2T 7 + 2T 8 , 

B 7 = T 7 — T s 1 

B 8 = — Ti + T 2 + T 5 - T 6 + 2T 7 + 2T 8 , 

C 2 = Vi - y 2 - y 3 , 

C 4 = -2Vi + \/ 3 + V 4 + 2V 5 , 

c 5 = v 4 - v 3 , 

C 6 = -Vi + v 2 + v 3 + v 4 + v 5 - v 6 , 

D 2 = -A 1 + A 2 - A 3 , 

D A = -2A 1 - A 3 - A A + 2A 5 , 

D 5 = A 3 - A 4 , 

D 6 = A 1 -A 2 + A 3 + A 4 -A 5 + A 6 , 
Ei — Si — S2 , 
H 1 = P 2 -P 1 . 

The expressions of the functions Pi(x) and p'^x) for the S° baryons can be obtained 
from the corresponding results of S + by making the replacement s -h- -u (s -h- d). 
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Figure captions 

Fig. (1) The dependence of the magnetic form factor Gm(Q 2 ) of the S + baryon on Q 2 
at s = 3.0 GeV 2 and M 2 = 2.0 GeV 2 , at several different fixed values of the arbitrary 
parameter j3. 

Fig. (2) The same as in Fig. (1), but for the electric charge form factor Ge(Q 2 )- 

Fig. (3) The same as in Fig. (1), but for the H~ baryon at s = 3.2 GeV 2 . 

Fig. (4) The same as in Fig. (2), but for the S~ baryon at s = 3.2 GeV 2 . 

Fig. (5) The same as in Fig. (1), but for the A baryon at sq = 2.6 GeV 2 and M 2 = 
1.8 GeV 2 . 

Fig. (6) The same as in Fig. (5), but for the electric charge form factor Ge(Q 2 )- 

Fig. (7) The dependence of the magnetic form factor Gm of the S + baryon on cos# 
at Q 2 = 1.0 GeV 2 and s = 3.0 GeV 2 , at several different fixed values of the Borel mass 
parameter M 2 . 

Fig. (8) The same as in Fig. (7), but for the electric charge form factor Ge- 
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